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We  shall  describe  a  procedure  which  calculates  the  Interior  bal¬ 
listic  properties  of  guns  by  numerically  integrating  the  partial  dif¬ 
ferential  equations  appropriate  to  the  dynamics  of  the  propellant  gas. 

The  equations  which  describe  the  Interior  ballistics  are  one- 
dimensional  along  the  gun  barrel,  assuming  a  uniform  cross  section. 

The  equation  of  motion  of  the  gas  in  Lagrangian  form  is 
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where  u  is  velocity,  t  is  time,  m  is  mass  per  unit  area,  P  is  pressure 
and  Q  is  the  artificial  viscosity  for  handling  shock  waves. ^  The  equa¬ 
tion  of  energy  conservation  is 
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where  E  is  Che  Internal  energy  per  unit  mass  of  the  gas,  V  -  1/p  is 
the  specific  volume  or  reciprocal  density,  and  B  is  the  rate  at  which 
energy  is  added  to  the  gas  by  the  burning  of  the  propellant.  The  equa¬ 
tion  of  motion  of  the  shell  as  it  moves  down  the  barrel  is 


du  , 
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where  is  the  mass  of  the  shell,  ug^  and  Pg^  are  the  velocity  and 
pressure  of  the  gas  at  the  shell  boundary,  and  A  is  the  barrel  cross- 
sectional  area. 

For  our  example  we  shall  use  a  large  5-inch/54  caliber  naval  gun 

2 

with  the  following  dimensions:  chamber  length,  Lo«100  cm;  cross- 
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sectional  area,  A-130  cm  ;  shell  mass,  M^-3.8  x  10  gm;  and  total 
shell  travel  distance,  L«600  cm.  The  propellant  used  is  90Z  nitrocel¬ 
lulose  and  10Z  other  volatiles,  and  it  is  in  the  form  of  small  cylin- 
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ders  with  diameter,  D-.69  cm;  solid  density,  pq-1.5  gm/cm  ;  and  total 

3 

weight,  VH9.0  x  10  gm.  The  propellant  contains  small  diameter  webbing 
which  we  shall  neglect. 

The  equations  of  state  of  the  gas  products  resulting  from  the 

burning  of  nitrocellulose,  mainly,  COj,  CO,  HjO,  H^,  and  N^,  have 
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been  determined  previously.  '  The  thermal  equation  of  state  is 

P(V-b)  -  nRT  ,  (4) 

3 

where  the  covolume  correction,  b».95  cm  /gm;  the  amount  of  gas  formed 
from  a  gram  of  nitrocellulose,  n- .04  gn-aole/gm;  and  the  gas  constant, 
R-8 3 . 15  bar*cm3/(gm-mole**K) .  The  caloric  equation  of  state  is 


E  -  P(V-b) /  (y-1)  , 


(5) 


where  the  ratio  of  specific  heats,  y-1. 25,  is  essentially  independent 
of  density.  Assuming  a  gas  ignition  temperature  (adiabatic  flame  tem¬ 
perature)  of  To"3000*K,  we  obtain  that  the  internal  energy  of  the  pro¬ 
pellant  is  Increased  by  about  Eo»800  cal/gm  upon  conversion  from  a 
4 

solid  to  a  gas. 

The  burning  law  for  a  cylindrical  propellant  of  diameter  fD,  where 
Olfll  as  the  propellant  burns  in  the  surrounding  gas,  is  empirically 


D  if  ■  ‘  BP  • 
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where  3  is  a  burn  rate  constant.  Neglecting  the  webbing  and  the  ends 

of  the  cylinders,  the  fraction  of  propellant  burnt  is 
2 

<(>-l-f  ,  and  it  follows  that 
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Thus  the  rate  at  which  energy  is  added  to  the  gas  is 


B  -  E„  if  -  E„  ?r  (1-*>1'’2 


(8) 


The  density  of  the  gas  as  a  function  of  4,  is 


p  -  W  $/[A(L+x)  -  (l-*)W/poJ  , 


(9) 


where  x  ■  x(4)  is  the  distance  which  the  shell  has  moved  down  the 


barrel. 


The  solution  is  obtained  by  numerical  integration  as  a  function 

of  position  and  time  of  the  finite  difference  equations  derived  from 

the  differential  equations  given  above. We  assume  i.he  shell  begins 

to  move  as  soon  as  it  overcomes  the  resistance  of  the  engraving  band 

4 

at  a  gas  pressure  of  F_,~400  bars .  .  This  determines  the  initial  gas 

conditions  at  which  we  start  our  calculation  as  p^*0.1  gm/cm  , 

Ej«350  cal/gm,  and  <p^».07.  Finally,  we  obtain  a  best  fit  solution  by 

specifying  the  burn  rate  as  B-.025  cm/(sec*bar) ,  which  happens  to  equal 

A 

the  experimental  value.  The  results  which  we  have  are  the  peak  breech 

pressure,  P  ,  the  shell  muzzle  velocity,  u  ,  and  the  shell  ejection  time, 
m  m 

t.  These  values  are  presented  in  Table  I  along  with  corresponding  ex- 
2 

perimental  values.  Also  a  history  of  breech  pressure  vs  time  is  shown 
in  Fig.  1. 

For  comparison  with  previous  work,  we  carry  out  the  same  ballistics 

calculations  using  the  representative  "Isothermal  solution"  described 
3 

by  Corner.  He  gives  expressions  for  P,  u,  and  x  as  a  function  of  the 
variable  f  and  the  "central  ballistic  parameter,"  M,  the  "form  factor," 

0,  the  "force  constant,"  X,  and  the  burn  rate,  B.  The  time  variable  is 
Introduced  through  Integration  of  the  time  dependent  equations  of  motion. 
We  obtain  a  bt'  fit,  but  not  unique,  solution  by  choslng  the  four  pa¬ 
rameters  to  give  the  correct  value  of  P  and  a  good  pressure-time  his- 
tory:  M-.93,  0**.6O,  A-240  cal/gm,  and  6*. 050  cm/ (sec* bar) .  Note  that 
3  is  not  in  agreement  with  the  experimental  value  above,  but  must  be 
used  for  consistency  with  the  other  parameters.  The  results  are  shown 
in  Table  I  and  Fig.  1. 

The  advantage  of  our  method  is  that  it  gives  the  details  of  the 
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gas  dynamics  along  the  gun  barrel  and  it  leads  directly  to  a  good  re¬ 
production  of  the  experimental  results  with  only  the  small  adjustment 
of  Che  burn  rate  parameter  and  it  is  capable  of  being  readily  modified 

to  account  for  small  corrections ,  such  as,  heat  losses  and  friction, 

6  3  4 
which  can  lead  to  improved  results.  The  earlier,  analytic  methods,  * 

such  as  the  one  listed  above,  are  often  detailed  but  are  usually  diffi¬ 
cult  to  modify  directly  and  usually  depend  for  accuracy  upon  the  some¬ 
what  arbitrary  adjustment  of  several  parameters. 


Table  I.  Sum*  17  of  Experimental  and 
Theoretical  Ballistics  Results. 
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